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Mixed-State Thermodynamics of Superconductors with Moderately Large 

Paramagnetic Effects 

Hiroto ADACHI * , Masanori ICHIOKA and Kazushige MACHID A 

Department of Physics, Okayama University, Okayama 700-8530 

Effects of Pauli paramagnetism on thermodynamic quantities in a vortex state, such as 
the specific heat C and magnetization M, axe studied using the quasiclassical Eilenberger 
formalism. We demonstrate that with an increase of paramagnetic depairing effect, the sigh 
of the curvature of the field dependence of C changes from negative to positive, and that 
the Maki parameter K2 becomes an increasing function of temperature. Our results provide 
a natural explanation for the unusual field dependence of C seen in CeCoIns in terms of the 
paramagnetic effect. 
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There has been much attention focused on heavy 
fcrmion superconductors, in particular on the newly 
found CeCoIns. 1 This material with a tetragonal sym- 
metry exhibits a first-order transition 2-5 at H C 2 both 
for H || c and H _L c. In a high field region for 
H 1 c, the possible existence of the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state 6-8 is highly debated cur- 
rently, because several experiments strongly indicate that 
the Pauli paramagnetic effect causes the first-order tran- 
sition, which may lead to the FFLO state. The pairing 



symmetry of this material on either d x . 



9, 10 



or gL,2_, 



is 

also controversial. One of the puzzling features concern- 
ing the pairing is that the specific heat C 9 ' u ' 12 obeys nei- 
ther the so-called Volovik effect 13 C oc \J~B expected for a 
line node case nor C oc B for a gapful case, 14-16 where B 
is the spatially averaged internal field. Instead CeCoIns 
shows a C oc i? 2 -like behavior. We point out that a sim- 
ilar behavior in C(B) is observed also in UBei3. 17 It is 
naively expected that the exponent 77 in C oc B v must be 
77 < 1, because each vortex carries a certain zero-energy 
density of states and C/C n > B/H c2 , where C n is the 
normal state specific heat. Thus, the observed C oc B 2 - 
like behavior is highly unusual and a puzzling issue to be 
solved. 

In a clean singlet type-II superconductor under a mag- 
netic field, there are two microscopic energy scales 18-20 
for the depairing effects: the Zeeman energy fiB for the 
paramagnetic depairing and the orbital depairing energy 
£b = vp/rs associated with the Doppler shift, where 
vf is a Fermi velocity, rs = y 4>o/2irB with <fio being 
the flux quantum. In materials with a small Fermi ve- 
locity such as heavy fermion superconductors, the role of 
the Zeeman energy compared to £3 cannot be neglected 
when discussing the thermodynamic properties. As for 
the phenomena near H C 2, the effects of Pauli paramag- 
netism on the mixed state have been studied by Adachi 
and Ikcda using the nonlocal Ginzburg-Landau (GL) the- 
ory. 5 To describe the thermodynamic quantities in lower 
fields we use the quasiclassical Eilenberger formalism. 21 

In this Letter, we clarify how Pauli paramagnetism 
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affects the mixed-state thermodynamics using the qua- 
siclassical Eilenberger formalism. Using the approximate 
method given in ref. 22, we demonstrate that in a su- 
perconductor with moderately large paramagnetic effects 
under a magnetic field, the specific heat varies roughly 
quadratically with the field, C oc B 2 , consistent with the 
experimental data in CeCoIns. 9 ' n ' 12 At a lower tem- 
perature with a stronger paramagnetic effect, our result 
shows a magnetization jump at H c2 , which is also ob- 
served in CeCoIn 5 . We give a simple thermodynamic 
argument to relate the positive curvature of C(B) to 
the magnetization behavior under the paramagnetic ef- 
fect, giving rise to the internal consistency between these 
quantities. Near H C 2, our numerical solution coincides, of 
course, with the results of the nonlocal GL theory given 
in ref. 5. 

First, let us briefly introduce our approach. We start 
with the Eilenberger equation 21 incorporating the para- 
magnetic effect 23 ' 24 (h = fee = 1); 

(2e M + iv n) /(e m , p, r) = 2g{e fl , p, r)u> p A(r), (1) 



where /(e M ,p,r), /t(e M ,p,r) 
0(£/*,P,r) 



-P, r)]*, and 



i/l — //t are the quasiclassical Green's 
functions. Here, v = vfP is a Fermi velocity, II = 
— iV + (27r/<f> )A is a gauge invariant gradient, T c is 
a transition temperature at a zero field. The fermionic 
Matsubara frequency e n — 2nT(n + 1/2) combined 
with the (renormalized) Zeeman energy is denoted as 
£ju = £n — i/^-Bj and the gap function is expressed as 
A p (r) = WpA(r) with the pairing function w p and the 
pair field A(r). In the following, the pairing state of a 
c?-wave (w p = 2^/2p x p y ) or s-wave (w p = 1) is assumed. 
Since the detailed shape of the Fermi surface does not 
change our main results essentially, we consider a quasi- 
two-dimensional isotropic Fermi surface. In the mixed 
state containing field-induced vortices, the pair field A 
can be expanded into each Landau level: 



A = A J2 d 

N=0 



(2) 
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where i> N = E~ = -oo O, "^? 1 e^H 1 /^, 

A = 1.764T C , C m = (i^/V^V^O 2 , ^ the N-th 
Hermite polynomial, and the lengths arc measured in 
units of Vb- The real constants Q and v are set to be 
( = 1/2 and v = (Sir 2 ) 1 ^ 4 , assuming a triangular vortex 
lattice. The difference in vortex configuration gives minor 
contributions. To solve eq. (1), we adopt the approximate 
method given in rcf. 22, which gives 

,OC Wmax 

/ = 2gw p / dp c- 2£ ^(A a N d N) • (3) 

N=0 

Here, the quantity is given by ajv = 

f~i H N (y+vm-Kc\) (\\\ 2 _ \ 2 ) /4^_ ( v + vm - \) 2 /2~wmx 

and A = (v y + iv x )p/rs- The above procedure uses an 
approximation similar to that used by Pesch, 25 in which 
the operator II does not pick up the spatial variation 
in g, but only that in A. For extreme type-II materials 
with a large GL parameter k (in CcCoIn 5 n ^> 10), 
this approximation provides a fairly good description 
for almost all experimentally relevant fields. Note that 
in contrast to the method by Pesch, our method can 
reproduce the first non-Gaussian term of the nonlocal 
GL free energy. 5 Recently, a simplified version of our 
method has been widely used 26 in slightly different 
contexts, in which the spatial degree of freedom is 
integrated out in advance. 

The expansion coefficients {djy} are determined by the 
following gap equation 



T/Tc =0.7 0.5 0.3 



T 



n>0 



2ttT 



~ r r 

^ = ^£^K(/ + (/*)*)>, 



n>0 



(4) 

where (• • • ) and • • • indicate the Fermi surface and spa- 
tial average, respectively. The strength of the param- 
agnetic effect in this Letter is measured using a para = 
M^orb/27rT c , where H OI \, = 0.561<I>o/27r£ 2 is the two- 
dimensional orbital limiting field and £ = vf/2itT c . 
If we define another parameter introduced by Maki as 
a M = V2H orb /H P (Hp = A /V2) , 18 ' 27 these two pa- 
rameters are related as % = 7.12a para . Throughout this 
Letter, we limit our discussion to the parameter region 
ctpara < 0.6 so that the possibility of either the vortex 
state formed by odd Landau levels 10 ' 24 ' 28 or the FFLO 
vortex state 10 is safely excluded, and we consider only 
the lowest Landau level. Furthermore, we define n as 
k- 2 = 47rAr(0)A 2 /F 2 rb , where N(0) is the density of 
states at the Fermi surface. 

To confirm whether our method can correctly capture 
the paramagnetic effect in the mixed state, we first study 
the magnetization. In the presence of the paramagnetic 
effect, the mixed-state magnetization is decomposed into 
the diamagnetic part Mdi a and the paramagnetic part 
M - 23 ' 24 

Jw para- 



4irM dh 



•">'.,.W .'/(/*("•, 



A) + f(w p A)*) 



B 



n>0 



1 + 9 



-2e n {l-g)) +c.c. 



(5) 




B/Horb 



Fig. 1. Field dependence of M<ji a in (i-wave superconductor with 
Opara = 0.6. Inset: corresponding data for a par a = 0. 



AttM, 



para 



4^M n (l+ (^)lm 



£<4 



(6) 



n>0 



where M n = 2tih 2 N(0)B is the Pauli param- 
agnetism in the normal state. In the clean limit 
without paramagnetism, the Maki parameter «2 = 
^[2p A d(AirM d{a )/dH]- 1 + 0.5| Hc2 with p A = 1.1596 is 
known to be a decreasing function of the temperature 
T. 29 ' 30 In the high-K case, we approximately have K2 oc 
[^(47rM c ji a )] B L / '^ , and the magnetization slope at H C 2 
gives a rough estimate of k 2 - The inset of Fig. 1 shows 
the field dependences of Mdi a in a d-wave superconductor 
without the paramagnetic effect (a para = 0). We can see 
that K2 is a slowly decreasing function of T in agreement 
with the results shown in refs. 29 and 30. Correspond- 
ing data with the paramagnetic effect (a para = 0.6) are 
shown in the main panel of Fig. 1. Contrary to the purely 
diamagnetic case, k 2 in this case is an increasing func- 
tion of T consistent with the observation in CeCoIns. 11, 12 
Our result should be compared with that of the dirty- 
limit analysis, 27 where K2 was shown to be an increasing 
function of T due to the paramagnetic effect. 

Next, we study the magnetization at a lower tempera- 
ture where a first-order transition at H c2 occurs. 5 ' 28 Fig- 
ure 2 shows the field dependences of Mdi a at T/T c = 0.2 
in a d-wave superconductor. In a low field region, the 
paramagnetic effect is irrelevant and all the data for a 
different a para value merge into a single curve. On the 
other hand, the slope dM d \ a /dB increases with increas- 
ing a para near iJ c2 , and for a para = 0.6, we find a dis- 
continuous jump in magnetization at H C 2 reflecting the 
first-order transition within the mean-field approxima- 
tion. Indeed, if we plot the free energy against the pair 
field strength near H C 2 (inset of Fig. 2), the free energy 
has a minimum at a nonzero value of (lAp) 1 ' 2 , implying 
a discontinuous jump in corresponding thermodynamic 
quantities. This explains the nearly discontinuous tran- 
sition at H C 2 observed in CeCoIn 5 . 2 ~ 4 We note here that 
we have also obtained essentially the same results also 
for the s-wave case. 

Now, we consider the specific heat C. To obtain C 
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numerically, we first calculate the entropy difference be- 
tween the normal and mixed states using 23 

V n>0 

| /t( mp A) + /(^ j ^ (?) 

where S n /V = 2tt 2 7V(0)T/3 is the entropy in the nor- 
mal state. Then wc perform a polynomial interpola- 
tion through the numerical data, and obtain C using 
C = TdS/dT. Figure 3(a) shows the field dependence 
of C/T at T/T c = 0.4 in a d-wave superconductor. In 
the absence of the paramagnetic effect (a para =0), the 
field dependence of C/T has a negative curvature, which 
coincides with the results shown in refs. 15 and 16. 
With increasing a para , however, the sign of the curvature 
changes from negative to positive. This is more evident 
in Fig. 3(b), where the same data are plotted against 
the normalized field B/H C 2(T). This field dependence is 
consistent with the observation in CeCoIns. 9 ' n ' 12 It is 
worth noting that this behavior can be seen irrespective 
of the structure of the pairing function. Actually, the 
same behavior is seen for the s-wave case, 31 as shown 
in the inset of Fig. 3(a). Also we would like to mention 
that although our result is based on an approximate so- 
lution and should be understood as a semi-quantitative 
one at most, an extensive numerical calculation 32 using 
the so-called explosion method supports the conclusion. 

To understand the microscopic origin of the behavior 
C(B) ex B 2 , it is convenient to relate the magnetization 
with C(B) through a thermodynamic Maxwell's relation, 

indicating that the T 2 -coefficient of the magnetization 
gives the slope of C(B)/T. Since the magnetization can 
be divided into the paramagnetic (M para ) and diamag- 



(Xpaia=0.6 0Cpara=0.4 0Cpara=0.2 Otpara=0 




0.2 0.4 0.6 0.8 1 



B/Horb 

Fig. 2. Field dependence of Mdi a at T/Tc = 0.2 in d-wave super- 
conductor for a para = 0, 0.2, 0.4, and 0.6. Inset: free energy near 
Hc2 vs pair field strength for a para = 0.6. The expression used 
for the free energy is the same as eq. (4) in ref. 23. 



netic (Mdia) parts, we can separately discuss the corre- 
sponding contribution to C(B) through the above rela- 
tionship. Figure 4 shows the temperature dependences 
of Mdi a and M para for a para = 0.6 at several fields. It 
is easily seen that M para changes like T 2 and this coef- 
ficient increases with increasing magnetic field, whereas 
M(ji a seems to be almost T-linear and has a small T 2 - 
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Fig. 3. (a) Field dependence of specific heat at T/T c = 0.4 for 
ci-wave case. Inset: corresponding data for s-wave case. Here, 
7n = 27r 2 A r (0)/3. (b) Data in (a) are plotted as a function of the 
normalized field B/H c2 (T). 
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Fig. 4. Temperature dependences of M para and M(jj a for 

C^p a r a — 

0.6 at several fields. Inset: field dependences of <9 2 M para /<9T 2 and 
<9 2 M dia /9T 2 at T/T c = 0.4 for a para = 0.6. 
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coefficient. This is clearer in the inset of Fig. 4, where 
the field dependences of d 2 M pal JdT 2 and d 2 M di JdT 2 
are shown; d 2 M paia /dT 2 is almost B-linear and gives 
the independence of (7(B), while d 2 M dia /dT 2 is almost 
field-independent and does not give a large field depen- 
dence. Thus, the main source of the behavior C(B) oc B 2 
is not the vortex contribution but the paramagnetic de- 
pairing of the Cooper pair that lowers H c2 consider- 
ably. It is interesting to point out that C(B)/T oc \f~B~ 
for a nodal gap case (ex B for a gapful case) at low 
temperatures indicates (3(B) oc \ j\f~B (oc const), where 
M(T) = M(0) + a(B)T + /3(B)T 2 + 0(T 3 ). Although 
there is no accurate measurement for M(T) to check it, 
the existing data of YNi 2 B 2 C 33 and MgB 2 34 seems to 
support it. It is desirable to confirm this internal consis- 
tency for other superconductors. 

We briefly discuss the data of Sr 2 Ru04, which is be- 
lieved to be a prime candidate for the triplet pairing. For 
H _L c, C(B)/T shows a strong positive curvature and 
changes smoothly into a \f~B~ behavior as H rotates to- 
wards the c-direction, 35 reminiscent of the change with 
apara shown in Fig. 3. This may be understandable if we 
assume either that the d-vector in Sr 2 RuC>4 is strongly 
locked within the basal plane 36 or singlet pairing symme- 
try. Upon rotating the field direction, the relative weight 
iJorb to the paramagnetic effect [iB, or a pa ra, changes, 
and C(B) /T should show positive to negative curvatures. 
Note that 37 (i) the in-plane magnetization curve exhibits 
a jump near H c2 at low temperatures, (ii) the Maki pa- 
rameter k 2 is an increasing function of temperature, and 
(iii) H c2 (± c) is strongly suppressed; these behaviors are 
similar to those of CeCoIn 5 . 

In conclusion, we have studied the mixed-state ther- 
modynamics of a superconductor with moderately large 
paramagnetic effects. Our numerical calculation based 
on the quasiclassical Eilenberger formalism revealed that 
with increasing paramagnetic effect, (i) the Maki param- 
eter k 2 (T) becomes indeed an increasing function of T 
and (ii) a C(B) oc B 2 behavior appears. It can be said 
that the field dependence of the specific heat and mag- 
netization is a good probe for measuring the strength 
of paramagnetic effects in type-II materials. Moreover, 
our demonstration provides a natural explanation for the 
C(B) oc B 2 behavior observed in CeCoIn 5 in terms of the 
paramagnetic effect. 

After this work was completed, we were informed 38 
that a C(B) oc B 2 behavior is a universal phenomenon 
for a system on the verge of exhibiting a first-order tran- 
sition of superconductivity; thus, it can also arise from a 
strong coupling effect such as antiferromagnetic fluctua- 
tions. 

We acknowledge fruitful discussions with T. Sakak- 
ibara, R. Ikcda, K. Tenya, T. Tayama, and K. Deguchi. 
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